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ABSTRACT 

The synchrotron emission observed from many astrophysical objects is commonly modelled 
to arise from relativistic electrons in an emission region occupied by a spatially homogeneous, 
but tangled magnetic field. However, we know that magnetic fields embedded in ionised gases 
tend to form flux ropes interspersed with regions of much lower magnetic field strengths. Here 
we develop a full description of the evolution of the energy distribution of relativistic electrons 
in a plasma divided into two distinct regions with different strengths of the magnetic field. 
Electrons are able to continuously leak from the low-field region into the high-field region. 
The model becomes fully analytic for physically reasonable assumptions. We show that such 
a leaky box model produces two distinct breaks in the electron energy distributions which 
give rise to three breaks in the resulting synchrotron spectrum. The spectral slopes in between 
the breaks are in general not constant and thus allow for significant curvature of the spectrum. 
These spectra are consistent with spatially resolved observations of the radio spectra of the 
lobes of radio galaxies. The exact form of the spectra depends on the adopted diffusion rate. 
The leaky box model significantly extends the time over which synchrotron emission can be 
detected at a given frequency compared to the usually assumed case of homogeneous magnetic 
fields. The spectral ages inferred for the electron population from standard techniques for the 
leaky box model are considerably younger than their real age. 

Key words: radiation mechanisms; non-thermal - radio continuum: general - methods: ana¬ 
lytical - galaxies: active 


1 INTRODUCTION 

Radio synchrotron radiation arises from relativistic electrons mov¬ 
ing in magnetic fields. The very nature of the synchrotron process 
enables us to infer the energy distribution of the electrons from th e 
observed radio spectrum (e.g. lKardashe\llT^|Pacholc^l 1 970l) . 
In many astrophysical objects this energy distribution is a power- 
law of the electron energy which arises from t he first order Fermi 
acceleration at shock fronts in gas flows (e.g. lB^ll978h . Over 
time, energy losses of the electrons due to the emission of syn¬ 
chrotron radiation as well as to the inverse Compton scattering of 
photons, mainly from the Cosmic Microwave Background (CMB) 
radiation, modify the shape of t 
ing a high-energy cut-off (e.g. 
high-energy cut-off gives rise to an exponential drop of the syn¬ 
chrotron spectrum at high frequencies. The position in frequency 
space of this drop is used to infer the ‘spectral age’ of the electron 
population, i.e. the time the electrons have spent in the magnetic 
field since they were accelerated to relativistic speeds. 

The standard spectral ageing technique assumes that the mag¬ 
netic field is homogeneous throughout the emitting volume and that 
its strength is constant for the spectral age of the electron popu- 
lation (e.g. iBur^ll977L IWinter et all Il98(i lAlexander & Leahvl 
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Il987t lAlexand^ Il987t IlTu et al.L Il992h . Adiabatic expansion of 
the emitting volume does not change the shape of the electron en¬ 
ergy distribution, but moves the high-energy cut-off to lower en¬ 
ergies as well as decreases the frequency of the associated spec¬ 
tral break. Consequently, the inferred spect ral ages are modified, if 
adiabatic expansion is taken into account iBlundell & Alexanden. 

s. 

Other effects that make spectral age estimates problematic 
are in-situ re-acceleration of the relativistic electrons and inhomo¬ 
geneous magnetic fields. Indications that thes e effects may pla y 
an important role were already reported by lAlexand^ il987h . 
The observation of radio spectra of spatially resolved lobes of 
radio galaxies with significantly more ‘curvature’ than expected 
from t he standard spect r al age i ng model emphasized these con¬ 
cerns iKatz-Stone et all 1 19931 : iRudnick et all Il994h . While re¬ 
acceleration clearly remains a possible explanation for the obser¬ 
vations, we concentrate in this paper on the effects of an inhomo¬ 
geneous magnetic fi eld. _ 

For Cygnus A iKatz-Stone & RudnickI H994I) found that the 
significantly curved radio spectra from any location within the ra¬ 
dio lobes, including the jets and hotspots, could be traced back to a 
single spectral shape by shifting the o bserved spectra in log (flux)- 
log (frequency) space. This result led iBlundell & Rawlings! i200(i) 
to develop a model for powerful radio galaxies in which a gradi¬ 
ent in the magnetic field strength along the lobes is responsible 
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for ‘illuminating’ different parts of a uniform electron population 
throughout the lobes. It is not clear how the specific gradient in 
the magnetic field required by the model can be maintained and 
the required conditions for the very fast diffu sion of the re lativis- 
tic electrons through the lobes may not apply iKaise 3 . I 2 OOOI 1 . Thus 
models with a magnetic field inhomogeneous on scales of the radio 
lobes may not explain the observations. 

On the other hand, it is well known that the strength of the 
magnetic field in the solar plas ma is far from homogeneous on 
small scales te.g. lZw^ ^ Il987h and the obser vations of filamen- 
tary structure in the lobes of radio galaxies (e.g. lCarilli et Iil ll99l[) 
may confirm the existence of simil ar fl ux ropes in these st ruc- 
tures. lWiita & Gonal-Krishn^ jl990l) and lSiah & Jl99d) de¬ 
velop a model with two separate regions of different magnetic field 
strengths which make up the synchrotron emitting plasma. The 
electron populations within the two regions also remain separate. 
They show that, when observed with a telescope unable to resolve 
the two regions, such a configuration can significantly alter the ra¬ 
dio spectrum compared to the case of a homogeneous magnetic 
field and hence complicates estimates of the spectral age. 

In the case where the strength of the magnetic field in the 
emitting volume follows a Gaussian distribution, the radio spec- 
trum does not devel op an exponential break in the GHz range 
jTribbld.[I^ll994 . Radiative energy losses of the electrons lead 
to a steepening of the radio spectrum, but since the loss history 
of individual electrons is very different, no clear break develops 
at low frequencies. The overall spectrum shows significant curva¬ 
ture ', consistent with the observations of iKatz-Stone et aP i 199.4 
and iRudnick et a DiHsi). 

The idea of a synchrotron emitting volume with t wo regions 
of diff erent magnetic fie l d stren gths was revisited bv lEilek et alJ 
<1994 and lAtovan et alJ <2Q0(]ll . In the first work the relativistic 
electrons are allowed to either move diffusively through the plasma 
or leak from one distinct region into the other. Both approaches 
lead to radio spectra with more than a single, exponential break. 
Thus they predict significant deviations of the estimated spectral 
ages from the real ages of the electron population. However, the 
mathematical complexity of these models makes them difficult to 
study their predictions in detail or to apply them to observational 
data. 

In this paper we develop a simplified version of the leaky 
box model which is based on original ideas of the transfer of 
cosmic rays from the Galactic disc to the G alactic halo (e.g. 
lKardash^ll962tlGinzburg & Svrovatskiil llj64[)/n^main differ- 
ence between our m odel and the work of lEilA et al 1 <1994 and 
lAtovan et alJ <2004 is that the rate at which electrons leak from 
the region with a low magnetic field strength into the region with a 
high magnetic field strength is independent of the electron energy. 
While this is a great simplification of the physics, it allows us to 
formulate an analytical solution for the evolution of the electron en¬ 
ergy distribution. We are then able to investigate the predictions of 
this simple leaky box model in some detail and find that it explains 
many of the observed properties of radio synchrotron spectra. 

Throughout the paper we assume that pitch angle scattering 
for the relativistic electrons is efficient and that a uniform pitch an¬ 
gle distribution is therefore maintained at all times. This behaviour 


' The predictions of the models presented in iTribbld h993t) and lTribbld 
il994t) for the curvature of the spectrum never appeared in print, but can be 
found at http://www.rfcgr.mrc.ac.uk/~ptribble/papers/preprints/color 


implies that o ur model is of JP- tvpe jjaffe & Perolat [197.4 . rather 
than KP-type <Kardashevill962tfPacholczvklll974 . 

In Sectionl^we derive expressions for the energy distribution 
of the electrons in both the low-field and high-field regions without 
making assumptions about the exact evolution of the volumes of or 
magnetic field strengths inside these regions as a function of time. 
In Sectionl^we make further simplifications and introduce an ex¬ 
plicit time-dependence to all relevant quantities which then allows 
the construction of a fully analytical solution. We then use this an¬ 
alytical solution in Section|3where we discuss the implications of 
the model for radio synchrotron spectra and for the standard spec¬ 
tral ageing estimates. We summarise our results in Section|^ 


2 DEVELOPMENT OF THE MODEL 

In our model the plasma containing the relativistic particles and the 
magnetic field is split into two separate, but spatial connected re¬ 
gions. The first region, region 1 or the ‘low-field’ region, of volume 
V\ (t) contains only a weak magnetic field with an energy density 
UQi{t). Correspondingly, the radiative energy losses of relativistic 
electrons in this region will be small. The strength of the magnetic 
field in the second region, region 2 or the ’high-field’ region, of vol¬ 
ume V 2 {t) is greater than that in region 1 so that MB2(f) > “Bl(f)- 
The radiative energy losses of the electrons in this region will be 
larger than those in region 1, but the synchrotron emissivity of re¬ 
gion 2 will also be greater. 

The volumes of regions 1 and 2 may depend on time. Thus the 
energy densities of the magnetic fields and the energy distribution 
of the relativistic electrons, even in the absence of radiative losses, 
will also be functions of time. For simplicity we assume that the 
ratio Vi (f) /V 2 (f) remains constant. It is straightforward to adapt our 
results to more complicated cases where Vi and V 2 have different 
temporal behaviour. 

In our model the relativistic electrons contained in region 1 
are allowed to leak into region 2. We do not consider leakage in 
the reverse direction. The main motivation for this assumption is 
the mathematical tractability of the problem. With leakage from re¬ 
gion 2 back into region 1, we cannot formulate an analytic solution 
to the relevant equations. Leakage in the reverse direction is also 
less interesting for practical reasons. Given that we expect most of 
the observable radiation to be produced in the high-field region and 
that the electrons will have a relatively short lifetime in this region, 
leakage from region 2 into region 1 will mainly influence the low 
energy end of the electron energy distribution in region 1. The rel¬ 
evant electrons will not contribute to the emission at frequencies 
we normally observe. Finally, the gyro-radius of a given electron 
will be smaller in region 2 compared to region 1. It will therefore 
be more difficult for electrons to diffuse from region 2 to region 1 
compared to diffusion in the opposite direction. 


2.1 Electron distribution in the low-field region 

The electron number density in region 1 is governed by the kinetic 
equation of form 

(Y,f)]-[3g(t)+^i (t)]«i(Y,f), (1) 

where 
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describes the dilution of the electrons by the expanding volume. We 
express the energy dependence of the electron density evolution in 
terms of the relativistic Lorentz factor, y. The function F (j,t) con¬ 
tains all processes which affect the energy of individual electrons. 

The leakage of electrons from region 1 into region 2 is de¬ 
scribed by the sink term hini. Note that the diffusion rate, hi , does 
not depend on y. In other words, of all electrons with a certain 
Lorentz factor y at time t, a fixed fraction diffuses during a time 
interval dt. The size of this fraction does not depend on y. In gen¬ 
eral hi will depend on the energy of the diffusing electrons. How¬ 
ever, we will show in the following that this restriction is necessary 
for an analytical solution of equation Q. Usually we will assume 
that the energy distribution of the relativistic electrons in region 
1 follows, at least initially, a power-law in y with a negative expo¬ 
nent. Thus there are many more electrons at low energies compared 
to high-energy electrons. Consequently, more low-energy electrons 
will diffuse into region 2. These low-energy electrons accumulate 
in region 2 and will eventually contribute significantly to the overall 
pressure inside this region while their contribution to the pressure 
in region 1 diminishes. At this point we would expect that further 
diffusion of low-energy electrons is suppressed, thereby introduc¬ 
ing an energy dependence in /ii. At this point our solution derived 
below becomes unphysical. However, again only the low-energy 
end of the electron distribution is affected which does not influence 
the synchrotron emissivity at common observing frequencies. So 
for practical purposes we can ignore this effect. 

We consider two forms of energy loss for individual relativis¬ 
tic electrons. The adiabatic expansion of Vi not only dilutes the 
electrons, but electrons also lose energy due to the work done dur- 
ing the expans ion. The adiabatic loss term is given by g (t) y (e.g. 
lLongaiilll994h . Radiative energy losses due to synchrotron radia¬ 
tion are proportional to y^, if the magnetic field is tangled on scales 
much smaller than the dimensions of the volume occupied by the 
plasma. The same proportionality holds for energy losses due to in¬ 
verse Compton scattering of CMB photons. We can therefore write 
for these combined loss processes fi (t)y^. All energy losses are 
then given by 

F{y,t) = fi{t)y^+g{t)y- (3) 


By noting that 

dni {y,i) _ dni (y,0 dni (y,r) dy 

dt dt dy dt ’ ^ ' 

we can now separate equation 0 into two new equations. The evo¬ 
lution of the Lorentz factor of individual electrons is described by 

^ = -fi{t)'i-g{t)y, (5) 

which has the well-known solution (e.g. lKardashe^ll962h 


yoG(LO) 

yQf{,G{t',0)fi{t')dF + l’ 

with 


G(f,ti) =exp 



( 6 ) 

(7) 


yo is the Lorentz factor at f = 0 of electrons with Lorentz factor y 
at time t. Note that the high-energy cut-off, yo.max^ for the initial 
distribution injected into region 1 at time f = 0, translates into a 
cut-off at 


Yl.max (r) 


YO.maxfz (t, 0) 

Yo,max/oG(T,0)/i(L)df' + l 


( 8 ) 
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at all later times t > 0. This high-energy cut-off is present even if 
Y0,max ^ 

The second equation resulting from the separation of equation 
Q describes the evolution of the electron density as a whole as 

= [2y/i (0 - 2g (0 - hi (f)] m (Y,f). (9) 

Because we restrict hi to be independent of y, we can solve this 
differential equation by separation of variables to give 

ni{y,t) = ni,oexp|^ [2y/i (?')-2g(f')-ill (f')] df'l 

77(Yl,max-Y), (10) 


where tii Q is the initial electron energy distribution injected into 
region 1 at f = 0. The Heaviside function H arises from the high- 
energy cut-off of the electron distribution, see equation From 
equation 0 "'6 nols that 

Y/iW = -^W-^, (11) 


and upon substituting for y/i, separating the integral in the expo¬ 
nential and using equation we find 


«i (yt) 


«i,o (y) G'^(t,0)exp 

77 (Yl.max ~ Y) ■ 



( 12 ) 


2.2 Electron distribution in the high-field region 

We will now derive an expression for the electron energy distri¬ 
bution in region 2, 112 (y,t). Consider electrons leaking from re¬ 
gion 1 into region 2 during a short time interval Ax at a time x 
with 0 < X < f. The energy distribution of these electrons in re¬ 
gion 1 is given by A«i (yx,x) = ni (yx,x)/zi (x)Ax. The Lorentz 
factor of the electrons at the time of leaking into region 2 is yx. 
Particle conservation requires Vi {i)hi (x) = V 2 (x )/!2 (x) and also 
Ani (yxjXjV'i (x) = An 2 (yx,"!) F 2 ('t)- Thus the energy distribution 
in region 2 of the electrons crossing from region 1 into region 2 at 
? = X is given by 

An2(Yx,'':) =«i (Yx,'i:)1i2(t:)Ax. (13) 

After arrival in region 2, the energy distribution of the elec¬ 
trons is governed by a kinetic equation similar to equation Q, but 
without a sink term, 

|■A«2(y.^) = ^{[/2(0 y^+«Wy] a«2(y,0 } 

-3g(r)A«2(y,0. (14) 

We use the same separation technique as in Section lrTI The evolu¬ 
tion of the Lorentz factor of individual electrons in region 2 is then 
given by 

yxG(t,x) 

y=—=-LL_L2_!-, ( 15 ) 

yx/'G(r',x)/ 2 (r')dF+l 

The high-energy cut-off for the electrons contributing to An 2 {y,t) 
can be found by substituting yx,max for yx in this expression. Clearly, 
yx,max is the high-energy cut-off of the electrons leaking from re¬ 
gion 1 into region 2 at time x and thus can be determined from 
equation by replacing t with x. 

The energy distribution of the electrons is governed by 

^An 2 (y, t) = [2y/2 (0 - 2g (f)] An 2 (y, t ), 


(16) 

















4 C.R. Kaiser 


which has the solution 

An2(Y,0=An2(Yx,^)(Y) G\t,x), (17) 

where we have used equation jl5L 

Using equation <13> . we can replace An 2 {yx,l) with an ex¬ 
pression containing the short time interval Ax. The complete en¬ 
ergy distribution of all electrons leaking from region 1 into region 
2 can then be calculated by replacing Ax with the differential dx and 
integrating, 

«2(YT)=/' G‘*(f,x)ni(Yx,x)i! 2 (x)dx. (18) 

7x„i„ V Y / 


The lower integration limit is governed by the evolving high-energy 
cut-off of the relativistic electrons in region 2 , 


Y 2 ,max (^5 


Yl,max G{t^X) 

Yl,max (t)G (f', X)/2 (f') df'+ r 


(19) 


Substituting for Yi.max ("i^) from equation and rearranging then 
gives 


Y 2 ,max (^ 7 "^) 


YO.maxG(r,0) 
YO.maxlr + 1 


( 20 ) 


where 

^=G(x,0) rG(r',x)/ 2 (f')dr'+ f G [t’,Q) f, [t') dt’. (21) 

Jo 

For a given Lorentz factor y we expect, at least at early times t, 
that Y < Y 2 ,max and in this case Xmjn = 0 in equation <18> . At later 
times radiative energy losses will have moved the high-energy cut¬ 
off, Y 2 ,max> to Lorentz factors below y. In this case, equation <20t 
is an implicit equation for the integration limit Xnjjt,. This integra¬ 
tion limit also removes the need for the Heaviside function in the 
definition of ni (Yxi't) in equation <12> . 

In addition to the electrons leaking from region 1 into region 
2 , some electrons may have been injected into region 2 directly 
at time t = 0. The energy distribution of these electrons will have 
the functional form of equation < 12 > with ni Q and hi replaced by 
^ 2,0 and 0 respectively. Also, the evolution of the Lorentz factor 
of individual electrons is given by equation jl5> with x = 0 and 
Yx = Yo- 


3 THE SIMPLE LEAKY BOX 

In this section we will simplify the general expressions for the elec¬ 
tron distribution functions derived above by suitable choices for a 
number of the parameters involved in the model. The aim is to build 
an analytical model which is simple to use, but allows us to study 
the generic properties of the synchrotron emission from a leaky 
box. 


3.1 Initial energy distributions and diffusion rates 

The exact functional form of the injected electron energy distribu¬ 
tions, no,l and no, 2 > depends on the mechanism responsible for the 
electron acceleration. Here, we will not consider distribution func¬ 
tions with a contribution from an initial injection of electrons into 
region 2. The synchrotron emission of such electrons is short-lived 
compared to that of electrons leaking from region 1 into region 2 . 
At times where Y 2 ,max (f 1 0) < Y for electrons with Lorentz factor y 
giving rise to synchrotron emission at observable frequencies, there 


is no contribution left from electrons initially injected into region 

2 . 

For electron acceleration at shock fronts, the most likely sce¬ 
nario for the lobes of radio galaxies of radio-loud qua sars, we ex¬ 
pect no I = wnYn ^ with 2 < p <3 and a constant no (e-g- lB3ll978t 
iHeavens & Drurvlll988h . By choosing p = 2 we. can significantly 
simplify the expressions for ni and n 2 and for the purpose of il¬ 
lustrating the basic properties of the leaky box model, we restrict 
our further discussion to this case. There are no analytic solutions 
of the various integrals in the expressions for n\ and n 2 for p ^2. 
However, the only effect of choosing /> > 2 is a general steepening 
of the spectral slopes discussed below. 

From equation <121 we then find 


«i (y,0 


noY ^G'*(t,0)exp 



H (Yl,max 


Y) • (22) 


The electron distribution in region 2 from equation <18> simplifies 
to 


«2 (y,0 





h 2 (x)exp 



dx. (23) 


The simplest diffusion rate we can adopt is hi = constant. We 
have already made the assumption that the time dependence of Vi 
is the same as that of V 2 - Thus we also have h 2 = constant from 
particle conservation. This choice then results in 

ni (y, t) = n^y-^G^ (L 0) exp (-/ri?) H (Yi,max - y) , (24) 


and 

«2 (Y.O = noy^^Cf{t,0) ^ [exp(-/;iXmin) -exp(-/iir)]. (25) 

hi 


3.2 Explicit time dependence 


The time dependence of all relevant quantities in the model is fixed 
by our choice for Vi it) and V 2 (t) as functions of time. We will 
concentrate on Vi (l) bearing in mind our earlier assumption that 
the behaviour of V 2 (f) as a function of t is the same. 

The evolution of Vi is governed by the fluid dynamics the 
plasma of the leaky box is subject to. In many practical applica¬ 
tions of the model we will have Vi <== f" with a>0 and so we set 

I/i(0 = Fi,o(^l + ^) , (26) 

where Vi q is the volume at f = 0 and Iq is a ‘scale height’ de¬ 
termining the speed of expansion of the volume. Substitution into 
equation ^ yields 


g(0 


a 1 

3 f-f ?o' 


while equation 0 gives 


G(f,q) 


ti+tQj 


(27) 


(28) 


The assumption of a magnetic field tangled on small spatial 
scales allows us to treat the field simply as a magnetic pressure 
with equivalent energy density u and adiabatic index F = 4/3. For 
our two regions we then have ui (l) == U 2 (l) (r). We set 

(f) = Mj o (1 + where the index j can be either 1 or 2 . 

ui Q and U 2 fl are constants and by definition ui q < j< 2 , 0 - Thus the 
radiative energy losses are described by 


( 29 ) 
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where aj is the Thomson cross-section, is the electron rest 
mass, c the speed of light and mcmb the energy density of the CMB. 
A useful result is the solution of the integral 


These settings fix the density of relativistic electrons in both 
regions. The synchrotron emissivity, i.e. the luminosity per unit vol¬ 
ume, of each region is then given by 



^[l-G'^(f,ri)]}, (30) 

with 

b = 5-3/fl 

d = 1-3/fl. (31) 


4 RADIO SPECTRA OF THE SIMPLE LEAKY BOX 


With the results from Sectionl^we can now study the properties of 
the radio spectra arising from the simple leaky box model. We com¬ 
pare these with the spectral properties of radio galaxies. We there¬ 
fore adopt model parameters appropriate for these objects. Note 
that the behaviour of the model may be quite different for a differ¬ 
ent choice of the model parameters. Here we only want to illustrate 
the general behaviour for the situation in radio galaxies. 

The model for the dynami cal evolution of the l obes 
of radio galaxies of typ e FRII ipanaroff & Rilevl 1 1974) by 
iKaiser & Alexand^ ^1997^ predicts that yolume elements of the 
lobe containing relatiyi stic electrons injec ted at the same time 
evolve with V with iKaiser et '^ Il997h 


4 + P 

''“re(5 + p)' 


(32) 


The model assumes a power-law distribution for the density of the 
gas the lobe is expanding into and — p is the exponent of this dis¬ 
tribution. In the following we use P = 1.5. The adiabatic index of 
the lobe material, Tc, is set to 4/3, since we are only interested in 
the magnetic field and the relativistic electrons, both of which have 
a relativistic equation of state. 

The expansion time of the lobes of radio galaxies is of or¬ 
der 10^ years and so we set fo equal to that time. The strength 
of the ma gnetic field in the lobes of Cygnus A is estimated as 
~ 100 uG fcarilli et al.Lll99lll and so we set 62 (f = 0) = 100 /jG 
for the strong-field region and fii (r = 0) = 62 (t = 0) /lO in the 
low-field region. The initial energy densities of the magnetic fields 
in the two regions is given by Uj q = B^{t = 0) / (Stt). We neglect 
any initial injection of relativistic electrons into region 2 and so the 
initial pressure in region 2 is po = “ 2 / 3 . For pressure equilibrium 
between the two regions we require p = (uj -fue) /3, from which 
we obtain Ue, the initial energy density of the relativistic electrons 
in region 1. We assumed earlier a power-law energy distribution for 
the elect rons injected into re gion 1 with an exponent p = 2. Thus 
we have iKaiser et a il ll997h 


no 


Me 

nieC^ 


In (yo.max) 


1 

-1 

YO.max 


(33) 


We set Yo,max =10®, but note that the exact value of this initial 
high-energy cut-off is not important. 

Finally, we choose the diffusion rate hi = 0.2/to and set 
h 2 = hi. The latter statement implies that regions 1 and 2 have 
the same volume. The energy density of the CMB is mcmb = 
4.3 X 10^*^ ergscm^^, equivalent to a negligible redshift. 


/ Yj,max 

ni{y,t)jv{j,t)iy, (34) 

where jv{y,t) is the emissivity of a single electron with Lorentz 
factor Y In our cal culations below we use the analytical expression 
for y'v presented in iGhisellini et alJ il988l) . 

The total monochromatic radio luminosity of regions 1 and 2 
combined is given by Py = [Vi + V 2 ) £tot,v = V'lEi.v 4- where 

£tot,v is the emissivity that would be inferred from observing fy 
under the assumption of a homogeneous emission region. It is in¬ 
teresting to note that 

V1+V2 V2 hi 

— r. -£tot,v = £l,v + 7r£2.v = £l.v + 7—£2,v- (35) 

Vi Vi *2 

From equations J24t and I2.5> we then see that only the normalisa¬ 
tion, but not the shape of the spectrum of the radio emission from 
the simple leaky box model depends on the ratio V 2 /V 1 = hi/h 2 . 
In other words, the spectral shape depends only on the diffusion 
rate between the two regions, but is independent of their respective 
volume filling factors. 


4.1 Electron energy distributions and spectra 

The radiative energy losses of the relativistic electrons in region 
1 cause the high-energy cut-off of the electron energy distribu¬ 
tion, Yl.max. to evolve according to equation Since the radiative 
losses of electrons in region 2 are greater than in region 1 and be¬ 
cause electrons constantly leak from region 1 into region 2, Yl.max 
is also the high-energy cut-off of the electron energy distribution 
in region 2. However, n 2 also contains a break at Y 2 ,max (pt)), de¬ 
fined by equation <20> . For y < Y 2 ,max (TO). radiative losses have 
not yet moved the high-energy cut-off of electrons, which leaked 
into region 2 at a time x = 0, below their current Lorentz factor. 
Thus the shape of the energy distribution is still that of the energy 
distribution initially injected into region 1. Also, note that Xjj,i„ = 0 
in this case. For Lorentz factors exceeding Y2,max(fi0)' 'tmin > 0 
and electrons which leaked into region 2 at times x < Xmin do not 
contribute anymore to the energy distribution. Due to this effect the 
shape of the energy distribution will be steeper between Y 2 ,max {p 0) 
and Yi,max than below Y 2 ,max (f, 0). 

Figure Q shows the evolution of the high-energy cut-off, 
Yl.maxi and of the break, Y 2 ,max(ti 0 )- The early evolution is very 
fast due to the strong energy dependence of the radiative losses. 
Later on the maximum Lorentz factors decrease only slowly. How¬ 
ever, this does not imply that the resulting radio spectrum remains 
constant, because the strength of the magnetic field and the density 
of the relativistic electrons continue to decrease. 

As discussed above, the high-energy cut-off, Yi.max^ gives rise 
to two breaks in the radio spectrum of the leaky box, one in the 
contribution from region 1 (hereafter referred to as Vi) and the 
second in the contribution from region 2 (hereafter V2). This is il¬ 
lustrated in Figure]^ which shows the radio spectrum at an early 
time, t = 0.5to = 5 x 10® years. At this point in time, about 10% 
of all electrons initially injected into region 1 have leaked into 
region 2. Due to severe radiative losses in region 2, the break 
frequency associated with Y2,max(L0) (hereafter V3) has already 
moved to very low frequencies. We can approximately predict the 
location of spectral breaks by using the fact that the synchrotron 
emission of relativistic electrons with Lorentz factor y moving in 
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Figure 1. Evolution of the maximum Lorentz factor in the two regions. 
Solid line: 72, max in region 2 (high-field region). Dashed line: yi,max in 
region 1 (low-field region). 



Figure 2. Radio spectrum of the simple leaky box at a time t — 0.5to. Long- 
dashed line: Contribution from region 1 (low-field region). Short-dashed 
line: Contribution from region 2 (high-field region). Solid line: Total spec¬ 
trum of both regions added together. The vertical lines show the locations 
of the expected spectral breaks. Long-dashed, vertical line: Spectral break 
(Vi) in region 1 due to Yi,max- Short-dashed, vertical lines: Spectral breaks 
in region 2 due to 72, max nt low frequencies (V3) and due to 7 i,max at 
high frequencies (V2). 

a magnetic field of strength B peaks around y^Vg/S, where the 
non-relativistic gyro-frequency is given by Vg = eBf (Inmpc) (e.g. 

iRvbickr&Lightinlf^ . 

The gyro-frequency of electrons in region 2 exceeds that in re¬ 
gion 1 by an order of magnitude. Thus, the overall spectrum at low 
frequencies is dominated by the contribution from region 2 where 
large numbers of low-energy electrons radiate at higher frequen¬ 
cies than their counterparts with the same Lorentz factor in region 
1. Beyond the break frequency V 3 radiative losses in region 2 have 
decimated the number of electrons, the spectrum of region 2 steep¬ 
ens and the flat-spectrum contribution of region 1 dominates the 
overall emission. Note that due to the high frequency tail of the 
single electron emissivity, jy, the spectra extend well beyond the 
break frequencies. This effect allows the contribution from region 
2 to dominate the overall spectrum at frequencies higher than V 2 . 

Figurel^shows the radio spectrum at t = to - Region 2 now con¬ 



Figure 3. Same as Figure P] but for time t = to. Note that the break at V3 
has now moved to frequencies below the lower limit shown here. 



Figure 4. Same as Figure|3 but for time t = 4to. 


tains 18% of all initially injected electrons. This fraction is enough 
to allow the emission from region 2 to exceed that from region 1 
at almost all frequencies. At frequencies between V3 (not visible in 
the Figure) and Vj the overall spectrum is steeper than the spectrum 
of region 1 alone due to the radiative losses in region 2. Flowever, 
at frequencies above Vi the overall spectrum is flatter than that of 
region 1. Also, the spectrum extends to higher frequencies, beyond 
the cut-off in the spectrum of region 1 . 

Finally, Figure^shows the radio spectrum of the simple leaky 
box at a time t = 4fo, when 55% of all electrons initially in region 
1 have leaked to region 2. The overall spectrum is now essentially 
identical to the spectrum emitted by the electrons in region 2. Al¬ 
though this spectrum is still steeper at frequencies below Vi than 
that of region 1 , it is clear that at frequencies above Vi, and in par¬ 
ticular at frequencies greater than V2, the emissivity of the leaky 
box is significantly greater than we would expect from electrons in 
a homogeneous magnetic field, i.e. from region 1 alone. 

4.2 Spectral indices 

The synchrotron spectrum arising from relativistic electrons with 
a power-law energy distribution within a homogeneous magnetic 
field can be described by a power-law with an exponential break 
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Figure 5. Colour-colour diagram for the simple leaky box with various dif¬ 
fusion rates h\. Solid line: Evolution of the spectrum for h\ — 0.2/ro, this 
is the case for which the spectra are shown in Figures |^to|4| Short-dashed 
line: h\ — 0.02/ro. Long-dashed line: h[ — 2/tQ. The dotted line shows the 
case of a pure power-law spectrum, while the dot-dashed line applies for a 
homogeneous magnetic field in the emitting volume. The poi nts are derived 
from observations of the hot spot s and lobes of Cyg nus A iRudnick et alJ 
11994 based on the obseiwations of ICarilli et a 1 I 1991 I) . 

at high frequencies where radiative losses have create a high- 
energy cut-off in the electron distribution. It has been pointed out 
iKatz-Stone et aliri993HRudnick et allll994l that the spectra mea¬ 
sured along the lobes of well-resolved radio galaxies show consid¬ 
erably more ‘curvature’ than predicted by this simple model. The 
spectral curvature is best illustrated in a radio ‘colour-colour’ plot. 

We dehne the spectral index between two frequencies Va and 
Vb with Va > Vb as 

a^=log(|)/log(^), (36) 

where Ea and Eb are the emissivities at the relevant frequencies. For 
observations at three frequencies, here we use 0.33, 5 and 15 GHz, 
we can calculate two spectral indices and then plot these against 
each other in a radio colour-colour diagram. 

In the case of a pure power-law energy distribution for the 
electrons without a high-energy cut-off or radiative losses we ex¬ 
pect 0 (q 33 = 3 . This is shown as the dotted line in the colour- 

colour plot in Figure|5| The track of a population of electrons evolv¬ 
ing in a homogeneous magnetic field is shown by the dot-dashed 
line. The evolution for this and any of the following tracks proceeds 
from the top right comer towards the left and down. The radiative 
losses move the high-energy cut-off to progressively lower energies 
and the resulting exponential drop in the spectmm causes the steep 
track in the colour-colour plot, since the spectral index between the 
higher frequencies always samples a steeper part of the exponential 
decay. 

Taking spectra a t various points along t he lo bes of several 
nearby radio galaxies. iKatz-Stone et alJ il993h and iRudnick et alJ 
i 19941 noted that these were mainly located in the colour-colour 
diagram between the dotted line and the dot-dashed line. The d ata 
points in Figurel^reproduce the results of iRudnick et allfTw^l for 
Cygnus A. Thus these spectra are inconsistent with spectral evolu¬ 
tion in a homogeneous magnetic field, even in t he case of inefficient 
pitch angle scattering teatz-Stone et allll993^ . 

Figurel^also shows the tracks for the simple leaky box model 
for various diffusion rates. Depending on the diffusion rate chosen. 


the tracks pass through the region occupied by the observational 
data in the colour-colour plot. The simple leaky box model is there¬ 
fore consistent with observations and can reproduce the observed 
curvature of the radio spectra of the lobes of radio galaxies. Note 
however the data points indicating the regions of the lobes with the 
flattest spectra in the upper right-hand comer of Figure|3 They can¬ 
not be explained with the simple leaky box model discussed here. 
They may require an energy distribution of the electrons flatter than 
Yq at injection time. 

The ‘plateaus’ seen in the evolutionary tracks through the 
colour-colour plot arise when one of the breaks or cut-offs in the 
energy distribution of the electrons passes by the lowest observing 
frequency. The first (upper) plateau is caused by the break at V 3 . Af¬ 
ter V3 has passed the lowest of the three observing frequency, the 
spectrum between them is increasingly dominated by the emission 
from region 2. The point where the contribution from region 2 starts 
exceeding that from region 1 moves from low to high frequencies. 
This effect therefore first leads to a steepening of the spectrum at 
the low observing frequencies (see Figure]^ and thus creates the 
plateau in the colour-colour plot. The first plateau occurs at flatter 
spectral indices for smaller diffusion rates. The plateau is also more 
pronounced for larger diffusion rates, so that it is difficult to see this 
first plateau for the lowest diffusion rate, hi = 0.02%, in Figure|3 
Both effects are due to the reduced importance of the emission from 
region 2 compared to that from region 1 for small diffusion rates. 
For hi = 0.2% the evolutionary track even crosses the line where 
“0 33 ~ 5 plateau, so that the spectrum is steeper 

at lower frequencies than at higher frequencies. This spectral be¬ 
haviour cannot be explained in the standard model of spectral evo¬ 
lution in a homogeneous magnetic field. 

The second (lower) plateau is associated with the break at Vi. 
This break affects the emission of region 1 only and so the plateau 
is more pronounced for cases where the contribution of region 1 to 
the overall emission is most important, i.e. for low diffusion rates. 
As expected, the evolutionary tracks for very low diffusion rates 
approach the line for evolution in a homogeneous field, because 
the emission of region 1 dominates. For high diffusion rates the 
evolutionary track approximates the line CICQ 33 = (Xj 3 , because the 
emission from region 2 dominates in this case and radiative losses 
at the high-energy end of the electron spectrum in this region are 
compensated by leakage from region 1. Also, radiative losses in 
region 1 do not significantly affect the spectrum in the range shown 
in Figurel^for this case, because the spectrum evolves rapidly along 
the plotted track. The position of the plateaus may be taken as rough 
indicators of time, since the break frequencies evolve in exactly the 
same way in each of the cases shown here. 


4.3 Spectral ageing 

In the standard model of spectral evolution in a homogeneous mag¬ 
netic field the position of the exponential cut-off in frequency space 
is used to determine the age of the emitting electron population. If 
the radio spectrum arises from a two-field plasma with leakage be¬ 
tween the two regions, then the results from spectral ageing may 
be misleading. In fact, we will show below that the standard spec¬ 
tral ageing approach may produce spectra that fit observations very 
well, but predict spectral ages that are considerably shorter than the 
evolution time of the two-field plasma. 

To study this effect we use the overall spectrum of the simple 
leaky box with hi = 0.2% and f = % shown in Figure|3 We then 
employ standard spectral ageing techniques to derive the spectral 
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Figure 6. Standard spectral ageing techniques applied to a synthetic radio 
spectrum of the simple leaky model. Solid line: Spectrum of the leaky box. 
Long-dashed line: Spectrum from standai'd spectral ageing technique with 
p — 2.3. Short-dashed line: Spectrum using a fixed exponent for the initial 
electron injection of p = 2.5. The large crosses show the five data points 
used in the fitting process. 


age of th is spectrum. For this analy sis we follow the meth od dis¬ 
cussed in lAlexander & LeahvHl98^ an d lAlexand^ ^1987^ . 

We assume that the spectrum produced by the simple leaky 
box model is observed at five radio frequencies. The flux at the 
two lowest frequencies, in our case 10 and 38 MHz, are used to de¬ 
rive the low-frequency spectral index, a. Since the spectrum aris¬ 
ing from the evolution of the electron distribution in a homoge¬ 
neous field predicts a power-law spectrum at low frequencies, the 
exponent of the underlying power-law electron energy distribution 
can be determined as p = 2a-|- 1. In our case we find p ~ 2.3. 
Using p and the ‘observed’ emissivity at 38 MHz we derive the 
strength of the magnetic field under the assu mption of mini mum 
energy conditions in the emitting plasma (e.g. lLongaiilll994) . We 
find Snjin = 22 /jG which is between the magnetic field strength in 
the high-field region, B 2 = 58pG, and that in the low-field region, 
B\ = 5.8/jG, at the time of observation, t = to. 

The strength of the magnetic field is usually assumed to be 
constant for estimates of spectral ages which simplifies the deter¬ 
mination of the high-energy cut-off in the electron energy distri¬ 
bution (see equationl^. This simplification also neglects the adia¬ 
batic energy losses of the electrons. With these assumptions we can 
now calculate spectra for different ages of the electron population 
and minimise their x^-deviation from the five ‘observational’ data 
points to find the best-fitting spectrum and hence the spectral age, 

^spec• 

Figure 1^ shows the result of the spectral ageing analysis ap¬ 
plied to the radio spectrum arising from the simple leaky box model 
as described above. Although the spectrum from the electron distri¬ 
bution evolving in a homogeneous magnetic field closely approxi¬ 
mates the spectrum of the leaky box, the spectral age derived is only 
3.8 X 10^ years compared to the real age of 10^ years. The standard 
spectral age determination technique underestimates the age of the 
relativistic plasma by a factor 2.6. Also, the exponent of the energy 
distribution of the injected electrons derived from the data points 
at low frequencies is too steep. This error can lead to a significant 
underestimate of the total energy content of the emitting plasma. 

In many cases no measurements of the radio flux of a source 
are available at low frequencies. In these cases the value of p is a 


guess. Often p = 2.5 is used. Figure|^also shows the best-fitting 
spectrum for this assumption. Again, the spectrum approximates 
the spectrum of the leaky box, but the derived spectral age is only 
1.7x10^ years, a factor 5.9 younger than the real age of the plasma. 

The reason for the underestimate of the spectral age is that 
the spectra of the leaky box model with p = 2 can evolve into a 
single power-law implying p > 2. In other words, the spectra of the 
leaky box model can mimic the spectra of electron populations with 
p > 2 evolving in a homogeneous magnetic field without a high- 
energy cut-off. Thus the standard spectral ageing technique arrives 
at an exponent p > 2 for the initial electron energy distribution and, 
since there is no sign of an exponential cut-off in the spectrum at 
high frequencies, the derived spectral ages are small. From Figure 
13 it is clear that this effect is particularly strong for high diffusion 
rates. 

A secondary effect is the move of the spectral cut-off to higher 
frequencies by the stronger magnetic field in region 2 as discussed 
in Section l4~n in connection with Figure |3 This effect also leads 
to an underestimate of the age of the electron population by the 
standard spectral ageing technique. 


5 SUMMARY 

We have developed an analytical model for the evolution of the 
electron energy distribution in a relativistic, magnetised plasma di¬ 
vided into two regions with different magnetic field strengths. This 
leaky box model is designed to mimic the evolution of a realis¬ 
tic magnetised plasma in which the magnetic field concentrates in 
magnetic flux ropes rather than being distributed homogeneously 
throughout the plasma. The relativistic electrons are allowed to leak 
from the region with a low magnetic field strength into the high- 
field region, but not back again. Also, we assume that the diffusion 
rate associated with this leakage is independent of the electron en¬ 
ergy. For simplicity we assume that the diffusion rate is constant 
and that the initially injected electron population has an energy 
distribution proportional to y In both regions the electrons are 
subject to energy losses due to adiabatic expansion, synchrotron 
radiation and inverse Compton scattering of the CMB. 

The evolution of the energy distribution of the electrons in the 
low-field region is identical to that of electrons in the commonly as¬ 
sumed homogeneous magnetic field, except for the integrated loss 
of electrons to the high-field region. As expected, the distribution 
develops a high-energy cut-off at Yi^max in terms of the Lorentz fac¬ 
tor. Because of the continued leaking of electrons to the high-field 
region and the faster energy losses there, Yl.max is also the high- 
energy cut-off in the high-energy region. However, the electron en¬ 
ergy distribution in the high-field region also develops a break at 
Y2,max < Yl,max- The break occurs at the maximum Lorentz fac¬ 
tor of those electrons that leaked into the high-field region at the 
earliest possible time which are therefore those electrons that have 
suffered the largest energy loss in the high-field region. 

The emission spectra of the leaky box model contain three 
spectral breaks associated with Yl.max and Y2,max' The spectral 
slopes in between the breaks depend strongly on the diffusion 
rate between the two regions and the time of observation. In gen¬ 
eral, the resulting radio spectra are more curved than those aris¬ 
ing from electron populations evolving in homogeneous magnetic 
fields. Therefore the spectra from the leaky box are consistent with 
observations of curved spe ctra in radio galaxies iKatz -Stone el's! 
ll99.3HRudnick et allll994i . More complex behaviour of the spec¬ 
tra is possible for the more complicated choice of a time-variable 
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diffusion rate and/or a different energy distribution for the initially 
injected electrons. 

We also applied standard spectral ageing techniques to the ra¬ 
dio spectra from the leaky hox model. We found that, while the 
standard spectral ageing model produces good fits to the spectra, it 
significantly underestimates the age of the electron population. The 
standard spectral ageing model derives the exponent of the assumed 
power-law energy distrihution of the radiating electrons from the 
spectral slope at low frequencies. The leaky hox spectra show sig¬ 
nificant steepening of the spectra even at these low frequencies and 
thus the standard models overestimate the value of the exponent. 
The discrepancies between the estimated spectral ages and the esti¬ 
mated exponent of the power-law energy distrihution and their real 
values is greatest for large diffusion rates between the two regions. 

Finally, the leaky hox model can considerably extend the time 
over which a given population of relativistic electrons can be de¬ 
tected at a given frequency. The energy of the electrons is pre¬ 
served in the low-field region. When they leak into the high-field 
region they light up and lose their energy quickly. However, more 
electrons of a given energy are available at later times to leak into 
the high-field region than for the case of a homogeneous mag¬ 
netic field. Thus synchrotron emitting structures remain detectable 
for longer than usually assumed. This may help to reconcile the 
long dynamic ages of buoyant radio plasma in galaxy clusters 
with the short rad iative ages derived using the standard spectral 
ageing model fe.g.lChurazov et allEoOlUMcNamara et^l200ll : 

iBruggen et alll2002li. 
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